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Reptiles world’wide are hit hard by the sixth mass extinc¬ 
tion (Bohm et al. 2013, Ceballos et al. 2017). While habi¬ 
tat destruction is the most important factor in biodiver¬ 
sity loss, the introduction of exotic species is a close second 
(Pimentel et al. 2001, Sax & Gaines 2003, Kraus 2015). 
When introduced species are related to natives, hybridi¬ 
sation may compromise the genetic integrity of natives, 
which may lead to unintended loss of population fitness 
(Fitzpatrick & Shaffer 2007, Allendorf et al. 2010). 
However, the fitness of local populations maybe increased 
if genes with an adaptive advantage are introduced, for ex¬ 
ample when local populations have become inbred or are 
below viable population size (Rieseberg et al. 1999, Ab¬ 
bott et al. 2013, Kelly & Phillips 2018). One well-moni¬ 
tored example of genetic rescue is found in Sweden, where 
an inbred snake population {Vipera berus) was rescued by 
introducing males from a genetically variable population 
nearby (Madsen et al. 1999, 2004). Hence, introduction of 
non-natives is sometimes conducted to boost native popu¬ 
lations, but negative effects such as loss of locally adaptive 
alleles and outbreeding depression should be carefully con¬ 
sidered when performing such genetic rescue (Weeks et 
al. 2011, Whiteley et al. 2015, Hamilton & Miller 2016, 
Russo et al. 2018). 

We illustrate the conservation complications related to 
anthropogenic hybridisation, using a case involving grass 
snakes (Natrix) in the Netherlands. Based on mitochondri¬ 
al DNA (mtDNA) and nuclear DNA, Kindler et al. (2017) 
distinguish western N. helvetica from eastern N. natrix and 
found that the two species meet in a narrow hybrid zone 
with very limited gene flow in Western Europe. The spe¬ 


cies occurring in the Netherlands is N. helvetica (Kindler 
et al. 2017, Stumpel & Ianssen 2017). North of the docu¬ 
mented hybrid zone, on the boundary between the Neth¬ 
erlands and Germany, the two species appear to be allopat- 
ric (Stumpel & Ianssen 2017). Sparse Natrix records, not 
identified to the species level, do occur in the region (de 
W ijER et al. 2009). 

Despite worldwide population declines in reptiles, 
Dutch reptile populations appear to be stable or growing 
(Ianssen & de Zeeuw 2017). The distribution of Natrix 
in the Netherlands has a long history of human interfer¬ 
ence, with intentional displacement by naturalists record¬ 
ed at least as far back as the second half of the 19* century 
(ScHLEGEL 1862), and may stretch as far back as the Neo¬ 
lithic (Lenders & Ianssen 2014). Recent records evidence 
that several populations in the west of the Netherlands 
(in the vicinity of populations 7 and 8, Fig. 1) were found¬ 
ed by (unauthorised) translocations of native individu¬ 
als from neighbouring areas (van der Lugt & Siebelink 
2003). Illegal release of exotic N. natrix derived from Italy 
and Turkey is thought to have occurred here as well (van 
DER Lugt & Siebelink 2003). Another introduction in the 
south of the Netherlands (in the vicinity of population 9) 
is thought to have concerned N. n. persa (Elzenga 1974, 
van Buggenum & Hermans 1986, van Buggenum 1992, 
Bugter et al. 2014). 

Morphological observations seem to corroborate the 
history of Natrix introduction. Snakes with two white- 
yellow dorsolateral stripes, a feature normally observed in 
N. n. persa (Kreiner 2007), have regularly been report¬ 
ed in the west of the Netherlands (vicinity of populations 
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7 and 8). The phenotype of N. helvetica includes lateral 
blocked stripes and is easy to distinguish from N. n. persa. 
A recent study from Switzerland reported the occurrence 
of mtDNA haplotypes of native {N. helvetica) and intro¬ 
duced (both N. helvetica and N. natrix persa) grass snakes, 
associated with the occurrence of offspring with diverse 
colour patterns (Dubey et al. 2017). 

Despite the difference in climate, it would not be sur¬ 
prising if exotic species would be able to establish viable 
populations in the Netherlands (Bugter et al. 2014). It 
is known that native grass snakes survive freezing con¬ 
ditions during hibernation in wetlands in both natural 
hibernacula as well as dry structures built by humans, 
such as rail embankments (van der Lugt & Siebelink 
2003). Furthermore, nesting opportunities for snakes 
are provided by manure heaps (Gunther & Volk, 1996, 
Lowenborg et al. 2010, 2012, Lenders & Janssen 2014). 
Although there is no indication of a population decline 
of Dutch grass snakes, grass snakes are indicated as vul¬ 
nerable on the Red List of Reptiles and protected by na¬ 
tional law (CBS 2005, VAN Delft et al. 2007, Janssen & 
DE Zeeuw 2017). Conservation efforts are taken to ensure 
survival of existing populations and facilitate dispersal 
among them. 



Figure 1. Sampling of grass snakes {Matrix) in the Netherlands. 
Large numbered dots are new localities and small dots are from 
Kindler et al. (2017). Blue reflects N. helvetica mtDNA haplo¬ 
types (clade E in Fig. 2), red (clade 4, for which the closest oc¬ 
currence is an introduced population in eastern Germany) and 
grey (clade 7, Cyprus and the west coast of Turkey, N. c. cypriaca) 
reflect M. natrix mtDNA haplotypes. 


Table 1. Sampling of nine Dutch Natrix populations. * Popula¬ 
tions where the presence of dorsolateral stripes was previously 
recorded. 


Pop. Location 

n 

Dorsolateral 
stripes (n) 

mtDNA 

N. helvetica / 
N. natrix 

1 

Fochteloerveen 

1 

0 

1/0 

2 

Marken 

3 

0 

3/0 

3 

Vaassen 

1 

0 

1/0 

4 

Asselsche Fleide 

1 

0 

1/0 

5 

Kootwijkerveen 

1 

0 

1/0 

6 

Flouten 

1 

0 

1/0 

7 

Alphen aan den Rijn* 

7 

6 

2/5 

8 

Krimpenerwaard* 

27 

10 

0/27 

9 

Brunssummerheide* 

1 

Unknown 

0/1 


Previous records of exotic Natrix individuals in the 
Netherlands consist of anecdotal and indirect evidence. 
To monitor the current genetic make-up of Natrix popula¬ 
tions in the Netherlands, genetic data may provide insight 
into where introductions have led to populations with ex¬ 
otic genotypes. The two grass snake species N. helvetica 
and N. natrix are represented by distinct mtDNA clades 
(Kindler et al. 2013, 2017, 2018). Therefore, mtDNA can be 
used to verify if multiple Natrix taxa co-exist in the Neth¬ 
erlands, if these snakes are introduced or native, and if con¬ 
servation efforts should consider the presence of multiple 
taxa and their potential hybridisation. Based on an mtD¬ 
NA phylogeography for Dutch Natrix, in combination with 
phenotypical information, we discuss natural and anthro¬ 
pogenic hybridization in the context of conservation. 

Tissue samples or saliva swabs from 43 Dutch snakes 
from 9 populations (Fig. 1) were collected with permis¬ 
sion from the conservation organisations Stichting Het 
Zuid-Hollands Landschap and Natuurmonumenten. The 
emphasis of sampling was on populations 7 and 8, which 
are suspected to include introduced snakes (Table 1). For 
each individual the presence of dorsolateral stripes was re¬ 
corded. DNA was extracted with the Qiagen DNeasy Blood 
& Tissue Kit. Fragments of two mitochondrial markers 
(ND4, 565 bp and cyt b, 899 bp) were amplified following 
Kindler et al. (2013). Sanger sequencing was done com¬ 
mercially at BaseClear B.V., Leiden, the Netherlands. Se¬ 
quences were edited in Sequencher 4.10.1 (Gene Codes 
Corporation, MI USA) and collapsed into haplotypes with 
FaBox (ViLLESEN 2007). Haplotypes were aligned against 
the NCBI website using BLAST (Altschul et al. 1990). 
Previously published haplotypes were identified and re¬ 
corded using haplotype numbers (Kindler et al. 2013, 
2017). As the sequences generated here are ~200 bp short¬ 
er, all possible matches with haplotypes (100% identity) are 
reported. One cyt b haplotype was entirely new, and was 
added to NCBI (MN585707). The new haplotype was com¬ 
bined with the datasets of Kindler et al. (2013, 2017) in 
Mesquite (Maddison & Maddison 2015). To determine to 
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which clade the new haplotype belongs, a maximum likeli¬ 
hood phylogenetic tree was built with RAxML 8.2.4, with 
1,000 rapid bootstrap replicates and a search for the best 
scoring maximum likelihood tree, treating each mtDNA 
marker as a separate partition and using the nucleotide 
substitution model GTR-i-G (Stamatakis 2014). 

Dorsolateral striping was present in 47% of individu¬ 
als from populations 7 and 8. No stripes were reported in 
any of the other populations studied, although dorsola¬ 
teral striping was reported previously in population 9 as 
well (Janssen 2009). The 43 Dutch snakes possess 3 ND4 
and 4 cyt b haplotypes. Ten individuals possess a N. hel- 
vetica haplotype (Fig. 2, belonging to clade E in Kindler 
et al. 2013, 2017) and 33 individuals a N. natrix haplotype 
(Fig. 2, clade 4 and 7, Supplementary Table Si). All haplo¬ 
types of N. helvetica were previously reported in western 
Europe (Kindler et al. 2013, 2017). When combining the 
information on the geographical origin of ND4 and cyt h 
haplotypes, the N. natrix haplotypes found in the Nether¬ 
lands occur naturally in southern and eastern Europe, up 
to central Germany; haplotypes previously found further 
west (in western Germany and Great Britain) are introduc¬ 
tions (Kindler et al. 2013,2017). One individual from pop¬ 
ulation 9 possesses a haplotype of N. n. cypriaca, which is 
endemic to Gyprus and the west coast of Turkey (Fig. 2, 
clade 7, Kindler et al. 2013). There is no correlation be¬ 
tween mtDNA type and the presence of dorsolateral stripes 
(Supplementary Table Si). 

We confirm the presence of three Natrix mtDNA clades 
in the Netherlands, typical of two recently recognized grass 
snake species N. natrix (clades 4 and 7; represented by the 


de3 


Clade 4 Austria, Bosnia-Herzegovina, Bulgaria, 
Croatia, Czech Republic, Germany, Greece, 
Hungary, Poland, Romania, Serbia, Slovakia, 
Slovenia, Sweden, United Kingdom 


- Clades 

- Clades 




Clade 7 Bulgaria, Cyprus, Greece, Turkey 


Clade B 
- Clade C 


- v.ictue u 


Clade E France, Germany, Swizerland, . 
The Netherlands, United Kingdom ^ 


Figure 2. Placement of Dutch Natrix haplotypes in the range¬ 
wide phylogeography of grass snakes by Kindler et al. (2013, 
2017). Dutch haplotypes were allocated to clade E {N. helvetica) 
and clades 4 and 7 {N. natrix). 


colours red and grey in Fig. 1) and N. helvetica (clade E; 
blue). Based on sampling throughout the Dutch Natrix 
range, only N. helvetica is thought to occur naturally in the 
Netherlands (Kindler et al. 2017). The origin of the N. na¬ 
trix introduction could not be connected to the previously 
suggested source in Italy. It appears that individuals from 
population 7 possess both native mtDNA from N. helveti¬ 
ca and exotic mtDNA from N. natrix, whilst individuals of 
population 8 only possess exotic mtDNA (Fig. i). Kindler 
et al. (2017) previously reported geographically restricted 
genetic admixture in the hybrid zone between N. natrix 
into N. helvetica, but noted the natural occurrence of asym¬ 
metric mtDNA introgression from N. natrix into N. helve¬ 
tica. This pattern of asymmetric introgression could reflect 
eastward hybrid zone movement, which would be in line 
with N. natrix ‘enclaves’ nested within N. helvetica territory 
(Gurrat et al. 2008, Toews & Brelsford 2012, Wielstra 
et al. 20i7a,b). However, the two N. natrix mtDNA clades 
found in the Netherlands occur far eastwards from the nat¬ 
ural hybrid zone and hence their presence in the Nether¬ 
lands cannot be explained by natural processes. 

No direct relation between exotic mtDNA and deviating 
phenotypes was found. However, the appearance of both 
striped and unstriped snakes in those populations possess¬ 
ing native and non-native mtDNA haplotypes suggests hy¬ 
bridisation occurs between native and non-native Natrix 
individuals. These findings are in line with previous find¬ 
ings in populations where both Natrix species are present, 
and offspring were observed to possess a wide range of 
phenotypes (Dubey et al. 2017). Unmonitored or uni¬ 
formed hybridisation between distinct species is not new 
and is especially problematic if species status is still under 
debate (e.g.; van Riemsdijk et al. 2017a, b). The complicat¬ 
ed issue of ‘genetic pollution, the replacement of native al¬ 
leles by non-native counterparts, should thus be taken into 
account in the conservation of Dutch Natrix (Meilink et 
al. 2015). ‘Genetic pollution should also be considered in 
the context of improving connectivity among Dutch Natrix 
populations: without species identification based on nucle¬ 
ar DNA, the spread of N. natrix alleles into the range of 
N. helvetica might inadvertently be promoted. To assess the 
impact of‘genetic pollution’ on Dutch Natrix we urge that a 
follow-up study incorporating nuclear DNA is conducted, 
focussing on the southwestern Dutch populations. 
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